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SO (3) and SO (2) decompositions of numerical channel flow turbulence are performed. 
The decompositions are used to probe, characterize, and quantify anisotropic structures 
in the flow. Close to the wall the anisotropic modes are dominant and reveal the flow 
structures. The SO(3) decomposition does not converge for large scales as expected. 
However, in the shear buffer layer it also does not converge for small scales, reflecting 
the lack of small scales isotropization in that part of the channel flow. 



1. Introduction 

Kolmogorov theory of fully developed turbulence has defined the main stream for al- 
most all theoretical and applied turbulence investigations since it appeared in the 1941 



flKolmogorov (194l| )) 

Kolmogorov made two important statements (i) there exists a tendency of any turbulent 
flow towards isotropization and homogenization of small scales fluctuations, (ii) there 
exists an inertial range of scales where the 'almost' isotropic and homogeneous turbulent 
fluctuations are characterized by a power law spectrum with a universal —5/3 slope. Both 

statements are connected and only partially correct. 

First, it is well established experimentally and numerically, Frisch (1995| ), that already 



in the ideal isotropic and homogeneous high-Reynolds numbers limit turbulent fluctua- 
tions cannot be characterized by only one single set of spectrum exponent, i.e. velocity 
fluctuations are strongly intermittent. Intermittency is the way to summarize the fact 
that the probability density of velocity increments, 5ru = (u{x + R) — u(x)) ■ R, over a 
distance R cannot be rescaled by using only one single scaling exponents for all R, see for 



example Frisch (1995 ). Second, in almost all relevant applied situation one is interested 
in those ranges of scale where turbulence statistics is neither homogeneous nor isotropic. 
In the following, as an example, we will discuss in detail the important case of channel 



flows. Recent experimental, parg fc War haft (1998 ), and numerical investigations, Pumir 



|(1996| ) ; pumir fc Shraiman (1995| ), have shown that the tendency towards the isotropiza- 
tion of small scale statistics of shear-flows is much slower than any dimensional prediction 
even at very large Reynolds numbers; in contrast to what is predicted by the Kolmogorov 
1941 theory, some observables like the skewness of velocity gradients, exhibit persistence 
of anisotropies. The two above issues are connected. One cannot focus on the issue of in- 
termittency in high-Reynolds number homogeneous and isotropic statistics without first 
having a systematic control on the possible slowly decaying anisotropic effects always 
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present in all numerical or experimental investigations. Similarly, the understanding of 
complex non-homogeneous and anisotropic flows cannot avoid the problem of intermit- 
tent isotropic and anisotropic fluctuations. 



Despite the many systematic theoretical attempts to attack intermittency in isotropic 
and homogeneous turbulence, the problem is mainly unsolved except for a class of toy- 
cases considering the passive advection by Gaussian and white in time velocity fields 
of scalar, IGawedzki & Kupiaincn (1995|); IChcrtkov et al. (1995), or vector quantities, 

Arad et al. (2000 ). Nevertheless, a lot of 



Gawcdzki fc Kupiaincn (1995 )j_ Chc rtkov et 
Vcrgassola (1996|) ; |Lanotte fe Mazzino (19961 ); 



phenomcnological and a nalytical progresse s have been don e by applying stand ard statis 



tical closure techniques, Kraichnan (1972 ); Lesieur (1987 ); L'vov et al. (1997 ), or more 



recent phenomenological tools borrowed from dynamical system theory like the fractal 





Benzi, Paladin, Parisi & Vulpiani (1984); 


Parisi & Frisch (1983 


); |Frisch (1995 


); 


Bohr, 


Jensen, Paladin & Vulpiani (1998]); 


Bcnzi, Biferale & Toschi (1998); 


Grossmann, Lohse 


& Reeh (2000). 



similar statistical attempts have been transposed to the understanding of 'non-ideal' tur- 
bulence, i.e., turbulence in all those situations when anisotropy and non-homogeneities 
play an important role in the turbulent production and dissipation. This paper is meant 
to partially fill the gap between the quantitative systematic methodology used in 'ideal' 
homogeneous and isotropic turbulence and the qualitative, ad-hoc, description in terms 
of 'structures' used in the 'non-ideal' wall bounded flows. We show how the decomposi- 
tion of statistically stable observable, like moments of velocity increments, in terms of 
the irreducible representations of the group of rotations in three dimensions, SO (3), and 
two dimensions, SO(2), allows a quantitative systematic characterization of the isotropic 
and anisotropic fluctuations. Moreover, a connection between the projections on differ- 
ent eigenvectors of the rotational groups and structures like 'hairpins' and 'streaks' is 
possible. 

Structures called 'streaks' have been thought to be the mai n signatures of wall bounded 



flows in the viscous sub-layer since the pioneering works of Kline et al. (1967 ) who ob- 
served the existence of extremely well organized motions made of region of low and high 
speed fluid, elong ated downstream and alternating in the span-wise direction. Later, in 



Kim et al. (197l|) , 'streaks' were reported to be the dynamical responsible of turbulent 
production in the viscous sub-layer. Similarly, 'hairpins' have been the main persistent 
structures o bserved experimenta lly, Head fc Bandyopanhyay (1981 ) ; Wallace (1982] ) , and 
numerically, Moin fc Kim (1985 ); Kim fc Moin (1986 ), outsid e the viscous layer, in the 



turbulent boundary layer. By mean of a conditional sampling, Kim fc Moin (1986 ) were 
able to show that these 'hairpin' shaped structures are associated with high Reynolds- 
shear stress and give a significant contribution to turbulent production in the logarithmic 
layer. 

More recently, [Toschi, Amati, Succi, Bcnzi fc Piva (1999 ); Bcnzi, Amati, Casciola, Toschi 



fc Piva (1999] ), started a first systematic investigation of the intermittent properties of 



velocity increments parallel to the wall as a function of the distance from the wall in a 
channel flow simulation. In this case a clear transition between the bulk physics and the 
wall physics was recognized in terms of two different set of intermittent exponents char- 
acterizing velocity fluctuations at the center and close to the channel walls. Still a firm 
quantitative understanding of how much these intermittent quantifiers can be connected 
to the presence of persistent structures is lacking. For instance, in Benzi et al. (1999| ) the 
different behavior of velocity fluctuations in the buffer layer was explained as a breaking 
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of the Kolmogorov refined hypothesis linking energy dissipation to inertial velocity fluc- 
tuations, i.e. an effect due to the different production and dissipation mechanism caused 
by the presence of strong shear effects close to the walls. Clearly, such a kind of issue 
can only be addressed by using systematic tools which are able to quantify the degree 
of anisotropy and coherency at difference scales and at different spatial locations in the 
flow. 

In this paper we propose to use the exact decompositions of the correlation functions 
in terms of the irreducible representations of the rotational group SO (3) (in the bulk of 
the flow) and in terms of the irreducible representation of the rotational group in two 
dimensions, SO (2), (close to the walls) in order to quantify in a systematic way the rel- 
ative and absolute degree of anisotropy of velocity fluctuations. Furthermore, we show 
how a careful analysis of the data allows also for a connection between some coefficients 
of the decompositions and the more common 'structures' observed by simple flow visual- 
ization. We show how the SO (3) decomposition, being connected to the exact invariance 
under rotations of the inertial and diffusive terms of the Navier-Stokes equations, is able, 
where applicable, to exactly disentangle universal scaling properties of the isotropic sec- 
tors from the more complex behavior in the anisotropic sectors. We also show how the 
SO (2) decompositions in planes parallel to the walls, is a useful analyzing tool in order 
to quantify the relative change of planar anisotropy by approaching the boundaries. 

The paper is organized as follows: In section 2 we review the main theoretical consid- 
eration about the importance of the SO (3) decomposition in the Navier-Stokes eqs. In 
section 3 we present a systematic analysis of the SO(3) decomposition in a numerical 
channel flow data base. We discuss the results with particular emphasis on the univer- 
sality issue, i.e. independence from the large scale effects, and on how one can use such 
a decomposition to quantify the relative importance of structures like 'hairpin' in the 
bulk of the flow. In section 4 we review the main findings which pushed us to apply the 
SO (2) decompositions in planes well inside the buffer layer, i.e. where the SO (3) decom- 
position cannot be applied due to the presence of the rigid walls, and we show how the 
SO (2) analysis allows us to clearly distinguish the existence of 'streaks' like structures in 
a statistical sense. Section 5 is left to comments and conclusions. 



2. SO (3) decomposition 

SO (3) - rotational invariance - is one of the basic symmetries of the Navier-Stokes 
equations. However, it is broken by the boundary conditions or by the driving force of 
the flow, both of which introduce anisotropy and also inhomogeneities. For the sake of 
completeness let us start, as an example, with the the SO(3) decomposition of the 2nd 
order most general velocity tensor depending only on one spatial increment R. 

C aP {x, R) = ((u a (x + R) - u a (x)(up(x + R) - u fj (x)) (2.1) 



It is easy to realize, Arad et al. (19995 ), that this observable can be decomposed in terms 
of the irreducible representations of the three dimensional rotational group which form a 
complete basis in the space of smooth second order tensors depending on one vector R: 

C a/3 (x, R) = Y, R)B q a P&)i ( 2 - 2 ) 

qjm 

The notation in ( |2.2| ) is borrowed from the quantum mechanical analogue, i.e. j — 0, 1, .... 
labels the eigenvalues of the modulus of the total angular momentum, L 2 ; m = —j, .., +j 
labels the eigenvalues of the projection of the total angular momentum on one direction, 
say y; q labels the different irreducible representations with a given j; and B^ m (R) are 
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the eigcnfunction of the rotational group in the space of second order smooth tensors. 
For example for the fully isotropic sector, j = 0, we have only m = and a simple 
calculation shows that there are only two independent irreducible representations in the 
isotropic sectors, i.e. the well known result, Monin & Yaglom (1975), that we need only 



two independent eigenfunction in order to describe any second order isotropic tensor. 
These two eigenfunction can be taken to be: 



and therefore the decomposition (|2.2|) in the isotropic sector assumes the familiar form: 



C a p{x, R) = a 1 . 00 (x, R)8 a .j3 + a2,oo(x, R)R a Rp 
In appendix A we list the comp lete set of B^f™ fo r the case of second order tensors. 



(2.3) 
For 

higher tensor ranks we refer to Arad et al. (19996). The main physical information is of 
course hidden in the dependence of the coefficients a q j m (x, R) on the spatial location, 
x, and on the analyzed scale, R. We aim at using the decomposition (2.2) as a filter able 
to exactly disentangle different anisotropic effects as a function of the spatial location 
and of the analyzed scale. In previous studies, the main interest was focused on the 
theoretical issue of the existence of scaling behavior for the coefficients a q j m (x, R) and 
on its possible dynamical explanation in terms of the 'foliation' of the Navier-Stokes eqs 
in different j sectors, Arad et al. (199S, 1999a, &|). The typical questions addressed were 
whether coefficients belonging to different j sectors have different scaling behavior (if 
any) and, in the case, which kind of dimensional estimate for scaling exponents in the 
anisotropic sectors one could propose. As for the issues of sc aling behavior, due to the 
limitation of small Reynolds numbers in the numerical case ( Arad et al. (1999a )), and 
to the limited amount of inform ation available on t h e tensorial structure of the velocity 
field in the experimental case ( Arad et al. (199S ); Kuricn et al. (200C )) only partial 



answers have been found. Among them, the most important is the strong universality 
shown by the isotropic sector as a function of the local degree of non-homogeneity (and 
anisotropy), i.e. the strong universality showed by the scaling properties of the coefficients 
a q, oo( x , R) as a function of x in non-homogeneous turbulence (Arad et al. (1999a)). On 
the other hand, in this paper we would like to also propagate the SO (3) decomposition as 
an appropriate tool to analyze, characterize, and quantify the non-universal large scale 
geometric properties of the turbulent flow. 



As an example we take numerical channel flow (Amati, Succi & Piva (1997); Toschi, 



Amati, Succi, Benzi & Piva (1999)) obtained by a lattice Boltzmann code running on a 
massively parallel machine. The spatial resolution of the simulation is 256 x 128 x 128 
grid points. Periodic boundary conditions were imposed along the stream- wise (x) and 
span-wise (z) directions, whereas no slip boundary conditions were applied at the top and 
at the bottom planes (y-direction). The Reynolds number at the center of the channel 
is about 3000. In our case of channel flow we assume that, due to the homogeneity in 
planes parallel to the walls, there is only a dependence on the height y of all statistical 
observable. The coefficients a q j m (x, R) carry two types of information: (i) Their scaling 



behavior 



(x,R) 



R C ^r, 



which at least for small scales and large Re is hoped 



to be universal, i.e., position and flow independent]!] and (ii) their absolute or relative 
magnitudes which clearly are non-universal, i.e., position x and flow type dependent. 
These ratios characterize what kind of structures the flow contains. These are time and 

f The issue of universality of sectors with j > is far from being trivial. A lack of universality 
may be due to the existence of infrared (IR) or ultraviolet (UV) dive rgences in the no n-local 
integral induced by the pressure terms in the Navier-Stokes equations, Arad et al. (1995). 
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ensemble averaged quantities, obeying the underlying Navier-Stokes SO (3) symmetry, 
and we consider them to be a more systematic tool for structures characterization than 
snapshots of vortex sheets, worms, swirls or contour plots of either the velocity or the 
vorticity fields. 

When analyzing higher order structure tensors C a /3,,, 7 (x, R) the decomposition of type 



(2.2) becomes cumbersome soon. Moreover, in most experiments the full tensorial infor- 
mation is not available anyhow. Therefore, one has to restrict oneself to an abbreviated 
form of the SO(3) decomposition of the velocity structure tensor, namely, the SO(3) de- 
composition of the longitudinal structure function. In this case, being the undecomposed 
observable a scalar under rotations, there exists only one irreducible representation for 
any j sector, i.e. the usual spherical harmonics basis set Yj m (R). We decompose the 
longitudinal structure function 



4 P) ' 



\x, R) = ( ((u (x + R) — u(x — R)) ■ rY) (2.4) 
as follows: 

S^(x,R)=J2S^,R)Y^{R). (2.5) 

jm 

We expect that when scaling behavior sets in (presumably at high enough Re) we should 
find: 

S^(x,R) ~a jm (x)R^ . (2.6) 

Again, the a^(x,R) carry both the scaling information a^ l (x,R) oc and their 

non-universal amplitudes. 



A practical problem with the decomposition (2.5) of (2.4) is that for x close to the 
boundaries the scale R is restricted to lengths smaller than the distance from the waljf]. 
More generally, R cannot exceed a typical distance over which non homogeneities are 



overwhelming. Therefore we will also perform a decomposition of (2.4) which obeys the 
weaker SO(2) symmetry, i.e. rotational invariance in a plane for fixed distance y from 
the wall, 

D^(y,R) = J £d^(y,R)exp(imct>). (2.7) 

m 

The orientation dependence in a plane reduces to the dependence on an angle cj>. Again, 
the d^} (y, R) carry both scaling and amplitude information. 

Let us notice at this point, that the SO(3) decomposition has its roots on the intimate 
structure of the Navier-Stokes eqs, i.e on the invariance under rotations of the inertial 
and dissipative terms and on the relative foliations on different sectors of the rotational 



group of the equation of motion of any correlation function Arad et al. (19996). On the 
other hand, do not exist closed equations for two-dimensional observable and therefore 
the SO(2) decomposition can only be seen as a powerful tool to exactly decompose 
any observable in a fixed plane as a function of isotropic and anisotropic structures 
in the plane itself. Clearly, such a kind of decomposition can teach us a lot in those 
regions, like at the border between the viscous and turbulent boundary layers, where 
strongly anisotropic but planar structures named 'streaks' are supposed to carry the 
most important dynamical information of the flow. 



t Because of the trivial remark that the analyzing sphere cannot touch the walls. 
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3. SO (3) analysis of a turbulent channel flow field 

In previous studies most of the attention was paid on the isotropic sector of the struc- 
ture function decomposition ( |2.5| ), i.e. on the behavio r of sjj^Jx, R) as a function of the 



center of the decomposition, x, and of the scale R. In Arad et al. (1999a ) it was showed 
that the isotropic projection enjoys much better scaling properties than the undecom- 
posed structure function and that these properties are robust with respect to the changing 
of the local degree of anisotropy, i.e. with respect to the center of the decomposition, x. 
These findings support the idea of universality of the isotropic scaling exponents. Very 
little was possible to say about scaling of the anisotropic sectors because of lack of spatial 
resolution; the only qualitative statement was that the scaling exponent of the j — 2 sec- 



tor was roughly 4/3, as predicted by the dimensional argument given by Lumley (1967) 



or by Grossmann et al. (1994) 



Here we want to concentrate also on the more applied question of how much the different 
projections, independently on their possible scaling properties, can teach us about the 
preferred geometrical structures present in the flow at changing the analyzing position 
in the channel. 

In figs. || and j§| we present the three different contributions we have in the j = 2 non- 
isotropic sector|| extracted at the center of the channel (y + = 160) and at one quarter 
(y + = 80) respectively. 

The relative size of the S 2 ,n(y + , R + ) for different m and fixed y + characterize the ge- 
ometry of the anisotropic structures on the corresponding scale R. E.g., for y + = 80 the 
(j = 2, to = 1) mode is very pronounced on smaller scales, see fig. ||. We associate this 
with the hairpin vortices and other structures which diagonally detach from the wall and 
which are projected out by Y 2 \. For a visualization of the Y 2m see figure [|. In the center 
the (j = 2, to — 1) mode is two orders of magnitude less pronounced than at y + = 80. 
Our interpretation is that the diagonal structures from above and below have equal and 
opposite contributions. 

The most pronounced structures in the center are those parallel to the flow direction, 
i.e., (j = 2, to = 2), see figure |[ 

Also at y + = 80 the structures parallel to the flow direction (mode (j = 2, to = 2)) are 
rather pronounced. At scales beyond R + rj 100 they overwhelm the diagonal contribu- 
tions (mode (j — 2, to = 1)). Therefore one is tempted to interpret R + ks 100 as the 
maximal size (in average) of the hairpin vortices. 

We re-did this type of analysis also for the S 2 ^l(y + , R + ) with very similar results. 
3.1. Higher order moments and the lack of isotropy at small scales 



The first question one may want to ask about the decomposition (2.5) is whether it con- 
verges with increasing j. We want to check this for an R in (stream- wise) flow direction, 
i.e., R = [9, (j>) — (7r/2, 0). As we can see from fig. ||, at small scales and in the channel 
center, where anisotropic contributions are small, the convergence is rather good. But 
away from the center (y + — 62) and in particular for large scales quality of the conver- 
gence become poor, see fig. [| Note that in any case the convergence is not monotonous 
as a function of the scale. This is a systematic quantitative way to understand the rate 
of isotropization toward small scales exhibited by this particular flow as a function of the 
distance from the wall. 

Another, even more informative way to quantify the rate of isotropization is to plot the 

(2) (2) 

ratio of each single amplitude S^{x, R) to the total structure function S L (x, R) with 

X The j = 1 sector is absent due to the symmetries of the structure functions chosen in this 
work. 
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R in the direction of the mean flow. In figs. ^ and || one can find the above quantities at 
the center of the channel y + = 160 and in the buffer layer y + = 62, respectively. What is 
very interesting to notice is that at large scales there are contributions from all resolved j 
sectors indicating as expected a lack of convergence of the decomposition at those scales 
and that in the buffer layer the relative ratio of the anisotropic sectors is much higher 
than what is seen in the center. Moreover, even more interesting, in the buffer layer, 
where due to the presence of a high shear one can imagine a statistically stable signature 
of anisotropic physics there appears a clear grouping of different sectors labeled by dif- 
ferent j indexes. Fig. ^| shows that projections with the same j but different m indexes 
have a qualitative similar behavior. Of course, these kind of comparison depends on the 
direction of the undecomposed structure functions (here taken parallel to the walls). 
Another possible test of the relative weights of anisotropies, free of the previous arbi- 
trariness, is to plot the ratio between the isotropic projection Sqq (x,R) and the other 
anisotropic projections for j > 0. Such a test is done in terms of quantities depending 
only on the separation magnitude \R\, and therefore measures the relative importance of 
anisotropies independently of the orientation. In Fig [5] we show, for example, the ratio 
between the sector (j, m) = (4, 4) and the isotropic sector (j, m) = (0, 0) at changing 
the analyzed height in the channel and for all R + . As it is possible to see, as expected, 
by approaching the wall (decreasing y + ) the ratio becomes larger and larger, showing 
clearly the importance of high j fluctuations in the sheared buffer layer. 
All the previous trends have also been found, amplified, by analyzing higher moments. 
For example, in figs. |l(] and |ll] we re-plot the same of figs. and ||| but for the fourth 
order structure functions. The fact that the previous trends are much more enhanced for 
higher order moments is a clear indication that anisotropy fluctuations are important but 
'rare', i.e. are connected to persistent intense fluctuations in a sea of isotropic turbulence. 



4. SO (2) analysis of a turbulent channel flow 

As extensively discussed in the previous sections, the SO (3) decomposition turned out 
to be extremely useful from both its theoretical background connected to the symmetry 
of the NS eqs and its ability to highlights statistical information as a function of their 
geometrical structures. On the other hand, the SO (3) decomposition suffers from some 
drawbacks when one wants to analyze the statistical turbulent behavior close to the fluid 
boundaries. This is due to the obvious fact that in order to perform the decomposition 
one needs to perform integrals over a given sphere, and therefore close to the boundaries 
the limitation of the sphere radius does not allow to extract any information but for a 
very limited (almost fully dissipative) range of scales. 

To overcome this problem we propose to use a decomposition in eigenfunction of the 
group of rotations in two dimensions, SO (2). The rational behind this idea is that the 
Navier Stokes equations obviously obey the SO(2) symmetry and for the channel flow 
also the geometry obeys this symmetry, once the rotation axis is chosen in the y direction. 
However, the mean flow breaks the SO(2) symmetry as it breaks the SO(3) symmetry. 
Nevertheless we will gain a tool being able to exactly decompose any two-dimensional 
observable in terms of fluctuations with a given property under two-dimensional rota- 
tions. In the region very close to the walls where very elongated 'streak' structures have 
been observed, the SO (2) analysis may help in understanding the relative importance of 
isotropic (in the plane) and anisotropic (in the plane) fluctuations. 

Another very important issue we want to address by using the SO (2) decomposition is 
connected to the recent findings by Foschi et al. (1999| ) of a different intermittent be- 
havior close to the walls (y + ~ 35) shown by longitudinal structure functions in the 
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stream-wise direction. These results were also connected to the breaking of the Refined 



Kolmogorov Similarity Hypothesis (RKSH) in the buffer layer, Bcnzi et al. (1999). Here, 
we show that by means of the SO (2) decomposition we are able to highlight the impor- 
tance of the streak like structures in determining this higher intermittent behavior. 
In this section we are interested only in observable in planes parallel to the walls and 
therefore the SO (2) decomposition of, say, the longitudinal structure function is defined 
as 

(y, R)=J2 d ™ (W. R ) ex P (M), (4.1) 

m 

where R is a two-dimensional vector lying in a plane at fixed y, and D < f\y,R) is the 
longitudinal structure function in the direction R. Due to the symmetry of the structure 



function we expect that only even ms will contribute to the sum in (4.1) 

In Fig. [l^ and Fig. |l^ we show the rate of convergence of the reconstructed structure 

function of order 2 as function of the maximum M contributing to the right hand side 



(RHS) of (4.1) and at two different distances from the wall, at the center (fig. |12|) and 
in the buffer layer (fig. |l3|). As it is clear, again we find a quite good convergence in the 
center. In the buffer layer, especially large scales are still far from being reconstructed 
even reaching M — 8. This is a clear tendency of formation of very large and intense 
structures in the strongly anisotropic buffer. These trends are even more pronounced for 
the fourth order moment as shown in Figs. ^ and ^|. 

In Fig. |r}] we show the absolute weight of different m-contributions for the second order 
structure function again in the center. It is interesting to notice that there is a clear mono- 
tonic organization of different contributions as a function of their isotropic/anisotropic 
properties, i.e. higher ms are less intense than lower ms systematically way at all scales. 
On the other hand, in the buffer layer, Fig. [H], there is a crossing of the m = 2 contribu- 
tion and of the m = 4 contribution at scales of the order of R + ~ 90. We interpret this 
crossing as the evidence of the formation of structures with typical size R + ~ 90 and 
with a preferred orientation given by the m = 2 eigenfunction. The m = 2 eigenfunction 
weights preferentially structures with a positive correlation in the stream-wise direction 
and negative correlated in the span-wise direction, i.e. exactly 'streak' like structures as 
those observed also in our numerical simulation by performing simple contour plots (see 
Fig. ^0j). As it is always the case, the above observed trends are even more intense for 
p = 4, 6... For p = 6 it even happens (not shown) that the dominant contribution at large 
scales is given by the m = 2 sector, proving, once more, the extreme departure from 
isotropy (in the plane) close to the walls. 

In order to quantify the departure from isotropy in each planes at changing the dis- 
tance from the wall we plot the ratios between the projections on the m = 2 sector 
and the isotropic sector (Fig. [l8]) at varying the distance from the wall and for some 
R + values. Fig. shows the same but for m = 4. It is interesting to notice, how there 
is a sharp transition for y + ~ 40 from an almost isotropic statistics {y + > 40) and a 
strongly anisotropic statistics (y + < 40), again the clear signature of the beginning of a 
"structured" buffer for (y + < 40) shows up. 

4.1. Review of near wall physics 
Let us now switch to the more statistically minded question connected to the existence of 



different intermittent properties close to the walls as previously reported in Toschi et 
K1999Q . This issue is connected to the general question whether in strong shear regions 

for a range of scales larger than the typical shear length, Lg = (e/S 3 ) 1 ^ 2 , one may have a 
different statistic transfer of energy than what expected at scales smaller than the typical 
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shear length. Different statistical energy transfer properties would have as a consequence 
also different intermittent properties of the velocity structure functions and, probably, 
the breaking of the RKSH if the shear lengths is small enough to be comparable with 
the dissipative length. Only in very high sheared region one can hope to have a very 
small Ls and therefore a sufficient range of scales with R > L$ where scaling laws can 
be investigated. Indeed in many different sheared flows, pnorato et al. (200C ) 
et al. (1996&); Gaudin et al. (1998| ), it has been found that in strongly sheared region 



intermittency increase a nd display univers ality (i.e. same exponents were measured in 
very different set-up). In Toschi et al. (2000 ) it was given a simple theoretical explanation 
of how intermittency should change in presence of strong shear. Considering the usual 
decomposition of the velocity field in its average and fluctuating part, v(x;t) = v(x) + 
v'(x;t), we get (from the Navier-Stokes equations) the usual Reynolds decomposition 



v'jdjV, 



D t v'i + Sij (x)v'j + v'jdjv'i 
Vjdj). The shear is defined as S. 



; = -d lP ' 



(4.2) 



djVi 



(x) and depends on the 



with D t = (dt 
mean flow geometry. 

In equation 4.2 the second and third terms of left hand side (LHS) will be of the 

1/2 

same order at a scale L$ ~ ( e /S 3 ) (shear length scale). For scale smaller than L$ 
it is the third term in eqn. [1.2] that will balance the energy dissipation, and hence the 
usual Refined Kolmogorov Similarity Hypothesis will hold: S p (r) ~ (e(r) p / 3 ) • r p / 3 . For 
scales larger than Ls it will be the second term in eqn. 4.2 that will balance the energy 
dissipation and hence: S p (r) ~ (e(r) p / 2 ). The validity of this second relation i n region 
of high shear values was already established by Benzi et al. (1999] ) . 
We want now to see how much one can say about this new 'intermittent' behavior close to 
the channel walls. In order to extract any quantitative information on scaling exponents 
in numerical simulation one needs to use the ESS technique, Benzi et al. (1993 , 1996a); 
Grossmann et al. (1997 ). ESS is based on the experimental and numerical observation 
that structure functions even at moderate Reynolds numbers show scaling in a generalized 
sense, i.e. by studying the relative scaling of one structure functions, say the second order 
structure function, versus any other. In particular, we want to verify and exploit that 
the following scaling holds: 

,( P)/C (2) 



D<*Hv,R)~(D$\v,R)) 



(4.3) 



where we have again limited ourself to the analysis of structure functions in the plane. 
In ( |4.3| ) we have explicitly taken into account the possi bility that the scalin g exponents 
depend on the distance from the walls. In particular, Toschi et al. (1999 ) showed, by 
analyzing the same data set, that there exist two distinguished set of exponents. One 
governing the scaling in the range of scales smaller than Ls (i.e. close to the center of 
the channel, in our case) which is given in terms of the usual isotropic and homogeneous 
set of exponents. The second governing the scaling in the sheared range of scales R > Ls 
(i.e. close to the walls in our channel simulation) which is given in terms of a much more 
intermittent set of exponents. 

In Fig. |2l] we show the ESS local slopes of the undecomposed structure function in the 
stream-wise direction in the center of the channel and the ESS local slope of the pro- 
jection on the m = sector always at the center of the channel for the moments p = 2 
versus p = 4. As it is evident, already the fully isotropic component (in the plane) is 
able to well reproduce the undecomposed observable and are both in good agreement 
with the isotropic and homogeneous scaling. Of course, the previous finding confirms the 
simple statement that at the center of the channel the shear length is formally infinite 
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and therefore the whole range of scales available is weakly affected by any shear effect. 
On the other hand in Fig. |2^ we show the same quantities of Fig. ^l] but in a plane well 
inside the buffer layer (y + = 37). As it is possible to see now the m — component 
does not reproduce the undecomposed observable, confirming the evident fact that here 
we are strongly anisotropic. Nevertheless, it is enough to add the m = 2 sector, i.e. to 



reconstruct up to M — 2 in the RHS of (4.1) to have a very good agreement with the 



more intermittent undecomposed structure functions local slope. This is a good quanti- 
tative evidence that as far as the new scaling properties are concerned the main effects 
is brought by these m — 2 'streak' like structures in the buffer layer. 
Still we need to understand the physics of these structures, why they are more intermit- 
tent, whether it is a coincidence or not that the new set of exponents coi ncides extremely 
well w ith what measured for passive scalar advected by a turbulent flow |Chavarria et al 



(1995 ). Nevertheless, we are confident that having a systematic way to analyze any 



isotropic/anisotropic two-dimensional/three-dimensional turbulent data set may help in 
further advance of the field. 



5. Conclusions 

A detailed investigation of anisotropies in channel flows in terms of the SO (3) and 
SO (2) decomposition of structure functions has been presented. Projections on the eigen- 
function of the two symmetry groups can be seen as a systematic expansions of structures 
as a functions of their scale and in terms of their local degree of anisotropy. 
We have used the SO (3) decomposition of structure functions at the center and at one 
quarter of the channel in order to have a quantitative tool to measure the relative im- 
portance of isotropic and anisotropic fluctuations at all scales. Close to the wall, the 
anisotropic fluctuations show strong effects induced from structure with the typical orien- 
tations of hairpin vortices. A partial lack of isotropization is still detected at the smallest 
resolved scales. 

The SO(2) decomposition in planes parallel to the walls allowed us to access also the 
viscous and buffer regions. In those regions, we have found that the strong enhancement 
of intermittency can be understood in terms of streak like structures and their signatures 
in some coefficients of the SO(2) decomposition. 

We think that the method presented here is beneficially applicable in all those cases 
where quantitative comparison and/or studies of anisotropic effects in different flows are 
needed (channel flows, boundary layers, homogeneous shear etc..) 

The application of similar decomposition to small scales observable like vorticity and 
energy dissipation would certainly be of great interest too. 

Having the possibility to control the anisotropic behavior is of great importance to im- 
prove LES of strongly anisotropic and inhomogeneous flows. 
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6. Appendix A 

In this appendix, we want to explicitly write down the SO (3) decomposition of the 
most general two-point velocity correlations in anisotropic turbulence. We consider the 
second order tensor involving velocities at two distinct points x x = x and x 2 = x + R: 

C a0 (R) = (-u a (R)u^(R + x)) . (6.1) 

where we have supposed that the statistics is homogeneous (but not isotropic) and there- 



fore the LHS of (6.1) depends only on R, the distance between the two points. Then, 
we can decompose C a>3 according to the irreducible representations of the SO (3) groups. 
Each irreducible representations will be composed by a set of functions labeled with the 
usual indices j = 0,1,.... and m = — j, corresponding to the total angular mo- 

mentum and to the projection of the total angular momentum on a arbitrary direction 
respectively. Moreover, a new 'quantum' index q which labels different irreducible repre- 
sentations will be necessary. It is easy to realize that there are only q = 1, 9 irreducible 
representations of the SO (3) groups on the space of two-indices tensor depending contin- 
uously from a three-dimensional vector, Arad et al. (19996 ). In particular, fixed j and m, 



the 9 basis tensor can be simply constructed starting from the scalar spherical harmonics 
Yj, m (X) plus successive application of the two isotropic operators R a and dp in order to 
saturate the correct number of tensorial indices. For example, the 9 linearly independent 
basis vectors which defines the irreducible representations in our case can be chosen as: 

Bf jm (R)^R^S a ^ 3m (R), 
Bf jm {R)^R-i +l e^d^ jm {R), 
B^ m {R) ee R-i[x a dP - R 11 d a ]<& ]m {R) , 
Bf jm (R) ee R-i- l e a ^R,<S> jm {R) , 
Bf jm {R)^R-^d a d^ 3m {R), 

B 6tjm (R) ee R-^le^R^d* + e a ^R^d^ jm (R) , 

Bf ]m (R) ee R~i{x a d p + R p d a )<S> 3m {R) , 

B^ m (R) ee R-i-^e^R^ + R?e a ^ R^d u ]<S> ]m (R) , 

B^ m (R)EER-^ x a R^ jm (R). 

Where for the sake of simplicity we have posed <F, m (R) = R^ Yj m (R) . As a results the 
most general second order tensor like (IO) can be decomposed as: 



9 



C«P(R) c ^m{R)B^ m {R) (6.2) 

j,m q=l 

where now the physics of the anisotropic statistical fluctuations must be analyzed in 
terms of the projections Cj jrriiq (R) in the different sectors. 
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Figure 1. Coordinate system in channel. Shown are stream- wise (a;), span- wise (2) and 

wall-normal (y) directions. 
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Figure 2. Log-log plot of s£l(y+, R + ) (*); S%l(y+,R+) (x) and S${y+,R + ) (+) as 
functions of R + at the center of the channel y + = 160. 
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functions of R + at y + = 80. 

t a 

Figure 4. Graphical representation of spherical harmonics (from left to right) |y 2,0 (#, <f>)\ , 

\Y 2 *\6A)\ and \Y 2 '\9,<t>)\. 
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Figure 5. Analysis of the convergence of the SO(3) decomposition: ratioj t M represents the 
ratio between the longitudinal structure function of order 2 in the stream-wise direction recon- 
structed up to (J, M) = (0, 0) (+); (J, M) = (2, 0) (x); (J, M) = (2, 2) (*); (J, M) = (4, 0) (□); 
(J, M) = (4, 2) (■) and (J, M) = (4, 4) (o) and the undecomposed structure function, at the 
center of the channel y + = 160. 
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Figure 6. Analysis of the SO(3) decomposition's convergence: ratioj t M represents the ratio be- 
tween the longitudinal structure function of order 2 in the stream-wise direction reconstructed 
up to (J,M) = (0,0) (+); (J,M) = (2,0) (x); (J, M) = (2,2) (*); (J, M) = (4,0) (□); 
(J, M) = (4, 2) (■) and (J, M) — (4, 4) (o) and the undecomposed structure function, at the 
center of the channel y + = 62 
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Figure 7. Ratio of each single (j, m) amplitude sj 2 ^ l (j/ + , J?) to the total structure function 

S}* (y + ,R) with R in the direction of the mean flow and y + — 160. The (j,m) indexes are: 
(0,0) (+); (2,0) (x); (2,2) (*); (4,0) (□); (4,2) (■) and (4,4) (o). 
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Figure 8. Ratio of each single (j,m) amplitude ^(y , R) to the total structure function 

S < i\y + ,R) with R in the direction of the mean flow and y + — 62. The (j,m) indexes are: 
(0,0) (+); (2,0) (x); (2,2) (*); (4,0) (□); (4,2) (■) and (4,4) (o). 
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Figure 9. Ratios between the sector (j,m) = (4,4) of the decomposition of S < i\y + , R + ) and 
the isotropic sector (j, m) = (0, 0) as functions of R, at changing the analyzed height in the 
channel: y+ = 160 (+); y+ = 125 (x); y+ = 92 (*); y+ = 80 (□); y+ = 62 (■); y+ = 48 (o) and 
y+ = 37 (.). 
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Figure 10. Ratio of each single (j,m) amplitude S < f} n {y + , R) to the total structure function 

S^i\y + ,R) with R in the direction of the mean flow and y + — 160. The (j,m) indexes are: 
(0,0) (+); (2,0) (x); (2,2) (*); (4,0) (□); (4,2) (■) and (4,4) (o). 
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Figure 11. Ratio of each single (j,m) amplitude Sj^ rl (y,R) to the total structure function 

S^(y + ,R) with R in the direction of the mean flow and y + — 62. The (j,m) indexes are: 
(0,0) (+); (2,0) (x); (2,2) (*); (4,0) (□); (4,2) (■) and (4,4) (o). 
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Figure 12. Analysis of the convergence of the SO(2) decomposition: ratiou represents the ratio 
between the longitudinal structure function of order 2 in the stream- wise direction reconstructed 
up to M = (+); M = 2 (x); M = 4 (*); M = 6 (□) and M = 8 (■) and the undecomposed 
structure function, at the center of the channel y + — 160 
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Figure 13. Analysis of the convergence of the SO(2) decomposition: rationi represents the ratio 
between the longitudinal structure function of order 2 in the stream- wise direction reconstructed 
up to M = (+); M = 2 (x); M = 4 (*); M = 6 (□) and M = 8 (■) and the undecomposed 
structure function, at y + — 37 



Probing structures in channel flow through S0(3) and S0(2) decomposition 21 




^* ********* K)to, „ 



_L 



_l_ 



_L 



_L 



_L 



_L 



_L 



0.5 

25 50 75 100 125 150 175 200 225 250 275 300 325 

R + 

Figure 14. Analysis of the convergence of the 50(2) decomposition: ratiou represents the ratio 
between the longitudinal structure function of order 4 in the stream- wise direction reconstructed 
up to M = (+); M = 2 (x); M = 4 (*); M = 6 (□) and M = 8 (■) and the undecomposed 
structure function, at the center of the channel y + — 160 
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Figure 15. Analysis of the convergence of the SO(2) decomposition: ratioM represents the ratio 
between the longitudinal structure function of order 4 in the stream- wise direction reconstructed 
up to M = (+); M = 2 (x); M = 4 (*); M = 6 (□) and M = 8 (■) and the undecomposed 
structure function, at y + — 37 
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Figure 16. Absolute weight ol different d„ (y + ,R + ) contributions for the second order structure 
function at the center of the channel y + — 160. The m's values of these components are: (+); 
2 (x); 4 (*) and 6 (□). 
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Figure 17. Absolute weight of different (y + , R + ) contributions for the second order structure 
function in the buffer layer y + = 37. The m's values of these components are: (+); 2 (x); 4 (*) 
and 6 (□). 
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Figure 18. Ratio between the projection on the m = 2 sector and on the isotropic sector 
m = as a function of y+, for R+ = 10 (+), 25 (x), 50 (*), 75 (□), 150 (■) and 250 (o). 
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Figure 19. Ratio between the projection on the m = 4 sector and on the isotropic sector 
m = as a function of y+, for R+ = 10 (+), 25 (x), 50 (*), 75 (□), 150 (■) and 250 (o). 
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Figure 20. Contour plot of the fluctuation of the stream- wise velocity at y + — 37. 
3 




Figure 21. ESS logarithmic local slopes of the undecomposed structure function in the 
stream- wise direction (+) and of the projection on the m = sector (x) as functions of the 
scale R + , for the moments p = 4 versus p = 2, at y + — 160. The dashed line represents the 
value 1.84 resulting from the experimental high-Reynolds numbers isotropic measurements 
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Figure 22. ESS logarithmic local slopes of the undecomposed structure function in the 
stream- wise direction (+); of the projection on the m = sector (x) and of the reconstruc- 



tion up to m„ 



= 2 



functions of the scale R + , for the moments p = 4 versus p = 2, at 



y — 37. The dotted-dashed line corresponds the best fit value, 1.52, for the ESS logarithmic 
local slopes of the undecomposed structure function in the stream-wise direction, the dotted 
line corresponds to the high-Reynolds number experimental isotropic value, 1.84, for the same 
quantity. 



